Searching for Dark Matter at the LHC with a Mono-Z 
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We investigate a mono-Z process as a potential dark matter search strategy at the LHC. In this 
channel a single Z boson recoils against missing transverse momentum, attributed to dark matter 
particles, Xi which escape the detector. This search strategy is related, and complementary to, 
monojet and monophoton searches. For illustrative purposes we consider the process qq — >■ XX^y 
where the Z boson is emitted from either the initial state quarks, or from the internal propagator. 
Among the signatures of this process will be a pair of muons with high pr that reconstruct to the 
invariant mass the Z, and large amounts of missing transverse energy. Being a purely elecroweak 
signal, QCD and other Standard Model backgrounds are relatively easily removed with modest 
selection cuts. We compare the signal to Standard Model backgrounds and demonstrate that, for 
large regions of parameter space, the signal would be clearly visible above background in future 
LHC data, allowing either new discovery potential or the possibility of supplementing information 
about the dark sector beyond that available from other observable channels. 
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I. INTRODUCTION 

There is now compelling cosmological evidence that 
the dark matter (DM) that appears to dominate galaxies 
and clusters of galaxies resides in the form of a gas of ex- 
otic weakly interacting elementary particles. While this 
has motivated direct searches for such candidates using 
underground detectors, and indirect searches attempting 
to find signatures of dark matter interactions in an as- 
trophysical setting, one of the most exciting possibilities 
involves producing and detecting such particles at the 
LHC. 

Although a multitude of dark matter candidates has 
been proposed in the literature (for some reviews see 
0]), by far the most popular class of model is that of 
a single massive particle species which interacts weakly 
with the standard model (WIMP). WIMP models pro- 
vide an appealing mechanism for dark matter produc- 
tion in the early Universe, compatible with the observed 
relic abundance, and, importantly, allow for the possibil- 
ity of direct and indirect detection through a variety of 
methods today. 

Assuming a WIMP-like dark matter candidate x, there 
are three types of x-Standard Model (SM) interactions 
that can be probed experimentally: annihilation (xx ~^ 
SM), scattering (x + SM ^ x + SM), and production 
{SM xx): which are of interest for indirect detection, 
direct detection, and collider searches, respectively. We 



consider here the last of these processes, specifically, the 
production of x at the LHC. 

If X couples directly to quarks, it will be produced 
through quark annihilation at the LHC, predominantly 
via the channel qq xX- Given the stability and weakly 
interacting nature of the DM, we expect x to leave the 
detectors unseen. Detection therefore requires a visible 
particle in the final state, against whose momentum some 
amount of missing transverse energy {^x) can be recon- 
structed. The dominant process which fits this criterion 
is dependent on the specific model Q- In many models, 
the simplest process that has both x and a visible par- 
ticle in the final state is the bremsstrahlung of a gauge 
boson during annihilation, qq — )■ xx+5ai*3e boson, where 
the gauge boson is emitted from either the initial state 
quarks or an internal propagator. Importantly, this pro- 
cess is generic and will therefore always be present, albeit 
with a rate which is model-dependent. 

The possibilities for gauge boson bremsstrahlung are 
the radiation of a single gluon (monojet), photon 
(monophoton), or electroweak gauge boson. The first 
two possibilities have been well studied. The hadronic 
monojets subsequent to gluon bremsstrahlung have been 
considered by several groups [8l-[l5| and are a promising 
dark matter search strategy. Similarly, monopliotons (7 
bremsstrahlung) can be used in dark matter searches, 
often in conjunction with jets [ll|, [3 [l3|, although con- 
straints are generally weaker than those derived from 



purely jet- based searches . 

The focus of this paper is to investigate colhder 
signatures of dark matter through purely electroweak 
bremsstrahlung. These processes can either lead to a 
a mono-VF |19j or mono-Z signal. Specifically, we con- 
sider a mono-Z signal, and highlight some unique kine- 
matical features of this channel, which make it an in- 
teresting and important complement to jet and photon 
based searches. Because of this unique kinematics, sig- 
nals may be distinguished from backgrounds even if rates 
are not as large as for other bremsstrahlung processes. As 
a result, depending upon the model, this new signal pro- 
vides either new discovery potential for dark matter at 
the LHC, or, equally interesting, information supplemen- 
tal to other observable channels to further pin down dark 
matter mo del- dependent parameters. 

We examine the expected signatures of the mono-Z 
process qq — XXZ at the LHC, relative to the SM back- 
grounds. We demonstrate these signatures by imple- 
menting a specific DM model in which the DM couples 
directly to quarks. This is used to demonstrate proof of 
principle for a mono-Z dark matter search, rather than 
being proposed as a fully self-contained particle physics 
model. However, many of the features discussed will be 
generically applicable in all WIMP models. 

Electroweak bremsstrahlung has recently received con- 
siderable attention in the context of dark matter annihi- 
lation and indirect detection [20. .33] . In certain mod- 
els, bremsstrahlung can play an important role in lift- 
ing a helicity suppression of the lowest order annihilation 
process, thus becoming the dominant annihilation mode. 
The possibility that lifting helicity suppression might en- 
hance electroweak bremsstrahlung associated with dark 
matter production at the LHC motivated our initial in- 
vestigations. However, we find that helicity unsuppres- 
sion negligibly affects rates in the kinematically accessi- 
ble detection regimes which we consider. Nevertheless, 
signals rise above standard model backgrounds. 

In Section |TT] we describe the dark matter - mono-Z 
signatures at the LHC, and outline the dominant SM 
backgrounds. In Section IIIII we introduce a simple DM 
model for which we examine the mono-Z signal, and cal- 
culate the production cross section. Here we also present 
a set of kinematic cuts designed to largely eliminate the 
SM backgrounds while preserving an observable number 
of signal events. Existing observational constraints on the 
model are discussed in Section llVi our main results are 
presented in Section |Vl and we conclude in Section |VT1 



II. LHC SIGNATURES AND BACKGROUNDS 

The XX^ production process is pictured schematically 
in Fig. [TJ In Section Hill we will introduce a specific DM 
model in order to explore this process in more detail. 
Here we will outline some of the general features of the 
XX^ final state, along with the relevant SM backgrounds. 
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FIG. 1. Generic electroweak bremsstrahlung process, qq — >■ 
XX^j which leads to a mono-Z signal at the LHC. 



A. The Zxx. Anal state 

Key to the discovery of the XX {Z ~> ff) final state 
from within the myriad of SM backgrounds is the correct 
reconstruction of a Z boson from the invariant mass of 
its decay products. We consider the muonic decay mode 
which, while having a low branching fraction (~ 3%), 
provides for a very clean invariant mass reconstruction. 
It also has the benefit of having few backgrounds relative 
to hadronic decay modes. 

The recoil of the Z against the heavy dark matter par- 
ticles results in decay muons with large transverse mo- 
mentum. The signal of interest will therefore be a pair of 
high pt muons, with an invariant mass that reconstructs 
to that of the Z, and a momentum sum which reveals a 
large amount of 

We simulate both the signal and background process in 
the MadEvent [s^l event generator. The MadEvent out- 
put is then fed into PYTHIA 35] to simulate higher order 
initial state and propagator radiation effects. A detailed 
detector simulation is beyond the scope of this work, as 
our goal is to demonstrate the potential for mono-Z dark 
matter processes to be observed above background at the 
LHC, rather than to calculate precise constraints on spe- 
cific dark matter models. 



B. Backgrounds 

The backgrounds for our process come from channels 
producing a dimuon pair and fx- The dominant back- 
grounds are the leptonic decays of gauge boson pro- 
ducing processes, specifically ZZ, ZW^, W+W' and 
tt ~> bbW+W. 

Production of Z + jets can also contribute a 
background through jet mismeasurement in the hadronic 
calorimeter. To test the importance of this background, 
we used MadGraph to simulate the dominant contribu- 
tion, Z plus a single jet. Hadronization was performed 
in PYTHIA, and fast detector simulation was carried out 
(for Z -t- jet background alone) in DELPHES [s^ (using 
ATLAS parameters), which simulates calorimeter smear- 
ing and reconstructs ^j^. As expected, this background 
was found to be sub-dominant to the other backgrounds 
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after the cuts described in Section llll Bl At leading order, 
and in the absence of full detector simulation and pile-up 
effects, the accuracy of our treatment of this background 
is obviously limited. Our simulations should, however, 
be accurate to within an order of magnitude of the true 
background, which is adequate for our purposes. 

As in the Z + jei case, mismeasurement of 6 jets will 
contribute to the spectrum for the process ti 
bbW~^W~ . This contribution is expected to be small 
compared to neutrinos from decay, and is therefore 
neglected. Similarly, mismeasurement of initial state ra- 
diation (ISR) in the form of gluon jets can contribue a 
background. Given the limited accuracy of PYTHIA in 
simulating these higher order process, we do not consider 
these effects. However, we expect their contributions to 
will be small, based on our simulations of the Z + jet 
background. 

The relative contributions of various backgrounds be- 
fore the implementation of the full set of cuts employed 
in this work can be seen in Fig. [2j Cuts on the invariant 
mass of the muon pairs ensure that NLO contibutions 
from 7 — fJ.'^ fi~ are negligible in these processes. The 
NNLO process 00 — ZZ can modify the ZZ background 
by up to 15% [37|; given the level of accuracy desired in 
this work, we neglect these contributions. These back- 
grounds are further reduced or eliminated through cuts 
described in Section llll Bl 

As evident from Fig.[2l the Z + jet background is sub- 
stantial. This same final state is also the dominant back- 
ground for the related monojet DM search channel. How- 
ever, we expect this background can be removed more 
easily for mono-Z's than for monojets. For a monojet 
search the invisible decays Z ^ vv provide a large fir 
background, with kinematics very similar to the XX + J^t 
signal searched for. In contrast, Z ~\- jets contributes to 
the mono-Z background through Z + jet /i"*"/^" + j^t, 
with arising only via jet mismeasurement. This is 
kinematically very diferent from our xx + ^ signal and, 
as we will show below, can be removed relatively easily 
with selection cuts. The sub-dominance of high cross sec- 
tion QCD backgrounds relative to electroweak processes 
is an appealing aspect of the mono-Z signal. 



III. THE MODEL AND EVENT SELECTION 

To illustrate the potential for observing a mono-Z dark 
matter signal at the LHC, we introduce a toy model in 
which this process has a significant rate. We will then 
detail event selection criteria that allow the backgrounds 
to be largely removed. 



A. An Example DM Model 

We take the DM to be a gauge-singlet Majorana 
fermion, Xi which couples to the quark doublet, Ql, via 
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FIG. 2. Modes contributing to full spectrum for signal 
and background for m-f, = 30GeV and m,, = 700GeV, after 
inclusive pr and invariant mass cuts. The signal cross section 
(— > XX Z) is calculated with a coupling that produces the 
observed relic abundance. 

the interaction term 

Ant = fudQL-n''x^ + h.c 

= fud {riuUi^ + riddL) Xr + h.c, (1) 

where fud is a coupling constant and t]'^ is a scalar field 
that that transforms as rj'^ ^ (3,2,1/3) under the SM 
gauge groups. (This model is a related to that of ref. [38j . 
modified such that the scalar is charged under SU(3)c-) 
Such couplings are also present in supersymmetric 
(SUSY) models, with x identified as a neutralino and 
rj a squark doublet. An obvious difference, however, is 
that we have no gluino analogue in our model. In some 
sense this model is analogous to a SUSY model in which 
the gluinos are too heavy to be kinematically accessible 
at the LHC. 

As a consequence, despite this model being substan- 
tially simpler than many SUSY models, both in cou- 
plings and free parameters, the LHC signatures presented 
in this work may still be of relevance for some SUSY 
searches (especially if the parameter space of more min- 
imal SUSY models is increasingly ruled out), perhaps 
providing a complementary signal to further constrain 
models. 

The interactions in Eq. [1] allow for direct annihilation 
of quarks into x P^-ir via t-channel and u-channel rj ex- 
change. Of interest to this work are processes to the 
next order in aw, in which a Z boson is radiated from 
the initial state quarks or the internal propagator. Con- 
tributing to the mono-Z process qq — XX'^ ^-re the three 
t- and u-channel diagrams shown in Fig. [31 

The Mathematica package FeynRules [3§| was used to 
formulate the Feynman rules for the model of Eq. [T] The 
rules were interfaced with the MadGraph package [s^l 
to calculate the scattering amplitudes for the processes 
in Fig. [3) These amplitudes were then input into the 
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FIG. 3. t-channel processes contributing to electroweak bremsstrahlung in annihilations to dark matter. Not shown are the 
three corresponding u-channel diagrams. 
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FIG. 4. Cross section for process pp — > XX^ CoM of 
14TeV as a function of DM mass. Red line (upper) corre- 
sponds to rrir) — 700GeV, blue (lower) to m,, = ITeV. Both 
cross sections calculated for f^d = 1, and using CTEQ6L1 
PDF's. 



MadEvent event generator, which calculated the LHC- 
relevant cross section pp XX^ for ^ given f^d, and for 
various choices of and . All cross section calcula- 
tions were performed in the proton center of momentum 
frame. The probability amplitudes were integrated over 
the CTEQ6L1 ^ parton distribution functions (PDF). 
Given that the LHC is a proton-proton collider, the va- 
lence quarks are all q distributions, while the sea-quarks 
are of course q and q distributions. 

Cross sections at 14TeV CoM are displayed in Fig.|3]as 
a function of m-^ , for values of relevant to electroweak- 
scale physics. 



B. Event Selection 

Now that we have a model which produces DM along 
with a Z-boson, we will examine how this channel may 
be detected at a hadron collider. While the backgrounds 
presented in IIIBI have rates much larger than our sig- 
nal, they can be substantially reduced with an educated 
set of cuts on measured events. We make cuts on the 
invariant mass of the muon pair within a 60GeV window 



centered on the Z mass, which greatly reduces the con- 
tribution from non-Z backgrounds, namely W^W~ and 
tt production. 

The presence of the heavy x in our signal process en- 
sures large amounts of This can be seen clearly 
in Fig. [21 which shows the number of expected collider 
events as a function of missing energy, in lOGeV bins. As 
expected, the number of signal events with large -^-p are 
at least comparable to all SM backgrounds. We choose 
a missing energy cut of $x > 150GeV to remove a large 
fraction of the background events, including the bulk of 
the Z + jet background. It is important to note that due 
to the very large cross sections for it and Z + jet before 
the implementation of cuts, the statistics in these two 
contributions lead to evident fluctuations at high in 

Fig.m 

The Z in the final state can be highly boosted by its 
recoil off the heavy DM particles; we therefore expect its 
decay products to have large pT- We apply the conserva- 
tive inclusive cut of pt > 50GeV on the muon transverse 
momentum (i.e. require at least one muon in final state 
with pT > 50GeV). 

A further consequence of the Z being produced rela- 
tivistically is that the muons from the decay process will 
be produced nearly co-linear with each other. This co- 
linearity ensures a low AR between the pair, where AR 
is defined to be 

AR = v/A02 + Ar;2, (2) 

where </> is the azimuthal angle and rj is the pseudo- 
rapidity of a particle in the detector. 

Figure [5] shows the ratio of signal cross section to the 
vP (Z — ^ fi'^ fj,~) background as a function of cut on max- 
imum AR, after pT and dimuon invariant mass cuts, for 
selective points in the model parameter space. This back- 
ground is useful for comparisons, as it is the dominant 
$T background in the region of interest. The signal to 
background is maximized for lower AR, with both cross- 
sections becoming equal around AR„iax ~ 1. To pre- 
serve signal events, we choose the conservative cut of 
AR < 1. This cut should effectively discriminate against 
the W~^W~ and ti backgrounds, which produce muon 
pairs with a broad AR distribution. The effectiveness of 
this choice of cut can be seen in the top right hand panel 
of Fig. [TTl which displays the same data as Fig. [5] but 
with full set of cuts applied, including the cut on AR. 
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FIG. 5. Ratio a^^z lozvv at 14TeV CoM and for ra^ — 
3 GeV, as a function of the cut on maximum Ai? = 
\J l^<fP + A772 between muons in the final state. Red line 
(lower) corresponds to = 700GeV, blue (upper) to rrir, = 
ITeV. 



The missing energy in the background Z-\-jet is a result 
of jet mismeasurement in the hadronic calorimeter, thus 
a large amount of fj, is present in events with high jet 
Pt- These types of events typically have a highly boosted 
Z, and subsequently produce muon pairs with lower Ai? 
separation than low $}rp events. 

We note that the detector simulation program 
DELPHES requires that muons counted individually be 
isolated within a cone of Ai? < 0.5. A consequence 
of these two factors is that in the majority oi Z + jet 
events with > ISOGeV, muons from Z decay do not 
pass DELPHES' isolation criterion and are subsequently 
rejected, reducing this background significantly. Aside 
from necessarily using DELPHES for producing the 
spectrum of the Z+jet background, we do not use detec- 
tor simulation in this work. Instead, we enforce Ai? > 0.5 
between the paired muons. 

Also, in this work we have been conservative in that we 
have not applied any jet related cuts, despite their obvi- 
ous utility in distinguishing real from jet mismeasur- 
ment. It is possible to further improve signal relative to 
QCD backgrounds by removing events in which jet and 
directions are correlated. Tagging of b jets may also 
be used to reduce the ti background. Below we show 
that such additional cuts, while no doubt useful, are not 
necessary to obtain an observable signal to background 
ratio. 
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FIG. 6. The minimum allowed coupling constant fud in or- 
der to satisfy the requirement that the contribution to the 
DM relic density from this model be less than or equal to the 
total DM relic density, flnAih^ — 0.11. Red line (lower) cor- 
responds to rrin = 700GeV, blue (upper) to m,, — ITeV. Note 
that the expansion parameter /^m remains perturbative for 
/ < 47r ~ 12.5. 



A. Freezeout 

We work in the context of the standard WIMP sce- 
nario, in which the DM was in thermal equilibrium in 
the early Universe up until the time of thermal freeze- 
out, at which point the relic abundance was set. For a 
given DM mass, we wish to choose values of the coupling 
constant fudi and 77 mass, such that the DM freezes out 
with the correct relic abundance. 

The process which kept x in equilibrium before ther- 
mal freezeout was qq — XX- The relic density of x was 
therefore determined by parameters fud and m^. Fol- 
lowing [41I, li^l , we used a semi-analytic solution to the 
co-moving Boltzmann equation, and the inferred value 
^DAih^ — 0.11 to place constraints on fud for given val- 
ues of m-^ and m^. Results are displayed in Fig. [51 If 
the coupling were any smaller than the constraints in 
Fig. [HI then the DM would have been overproduced in 
the early universe, yielding an abundance greater than 
that observed today. On the other hand, if it were any 
larger, then the relic abundance would be smaller than 
observed. If there are other DM species contributing to 
the relic abundance, then the constraints on fud serve as 
lower limits, since the DM candidate under consideration 
need not contribute the entire relic abundance. 



IV. MODEL CONSTRAINTS 

When calculating LHC signals, we adopt model pa- 
rameters that produce the correct DM relic abundance. 
We also ensure the parameters are in accordance with 
current direct detection and collider bounds. These con- 
straints are described in detail here. 



B. Direct Detection 

In the model under discussion, quarks couple to the 
beyond-SM sector via a qxr] vertex with strength fud- 
Consequently, care is required to avoid direct detection 
constraints. The operator in Eq. [T] allows for x-quark 
scattering via the s and w-channel 77 exchange graphs in 
Fig. [3 which can in turn be related to x-nucleon scatter- 
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FIG. 7. Fundamental processes contributing to DM-nucleon 
scattering. 
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FIG. 8. The spin-independent (SI) x— nucleon scattering cross 
section. Red line (upper) corresponds to m,, = 700GeV, blue 
(lower) to rrir, = ITeV. Shown in dashed is the upper limit on 
the SI cross section from XENONIOO experiment. 



mg. 

The couplings in Eq. [1] Fierz transform to couplings 
to nucleons that have both spin-dependent and spin- 
independent contributions. The strongest constraints on 
our model come from the spin-independent limits from 
the XENONIOO experiment 43.], which looks for excita- 
tions in ultra-cold liquid Xe resulting from DM scatter- 
ing events. We performed the calculation of the SI X" 
nucleon cross section in the current model using the Mi- 
crOMEGAs [i^l software package, taking the Lagrangian 
in Eq. [1] as input. The calculation was done for values 
of fud that produce the correct relic abundance (Fig. [B]) 
for various values of and a range of DM masses. The 
results are displayed in Fig. [8l alongside the upper bound 
on the cross section allowed by XENONIOO. Clearly the 
model parameters considered in this work are allowed by 
the XENONIOO constraint. Note that if there are addi- 
tional DM particles, and the DM candidate considered 
here is not required to contribute the entire relic density, 
the curves in Fig. |5] denote lower limits on the scattering 
cross section. 



C. Collider Constraints 

Through the early part of 2012, roughly 5 fb~^ of data 
have been analyzed by both the ATLAS and CMS col- 
laborations. Thus far no significant discrepancies with 



Standard Model predictions have been found. Discrep- 
ancies not seen at present include large missing energy 
signals and new particles such as those predicted in nu- 
merous extensions of the SM. 

The absence of novel signals has put ever tighter 
constraints on models of physics beyond the SM. As 
discussed in Section IIII Al the model described in the 
present work is purely phenomenological, but does pos- 
sess some similarities to a SUSY model (with obvious 
differences being an absence of giuinos and only first gen- 
eration quark couplings). This makes it somewhat diffi- 
cult to directly compare existing bounds with our model. 
However, we have chosen values of parameters such as the 
rj masses and rj — x niass differences which are not cur- 
rently ruled out by squark mass and squark-neutralino 
mass difference constraints, respectively, for a simplified 
model spectra (SMS) of SUSY. 

Examples of these constraints for an SMS from CMS 
data are given in [45 - 48, ]. The SMS are motivated by 
popular SUSY frameworks such as the constrained min- 
imal supersymmetric standard model (cMSSM) and the 
general gauge mediation model (GGM). Of particular rel- 
evance are the constraints derived in [461] . in which the 
authors consider a simplified model with a heavy gluino 
and two squark generations. This very closely resem- 
bles the model at hand, the only effective difference be- 
ing the number of scalar degrees of freedom present in 
the model. A lower bound of 780GeV was placed on the 
squark mass, which maps to the constraint > GOOGeV 
in the model considered here. Conservatively, we only 
consider r] masses above 700GeV. 



V. RESULTS 

The results of our simulations can be seen in Figs. |9]- 
[T3l These figures sample a range of model parameters, 
LHG energies and integrated luminosities, and use values 
of fud adherent to the constraints in Fig. [6l Plotted are 
expected signal and background events per lOGeV bin 
as a function of missing energy and the pt of the /i^, 
after the application of cuts outlined in Section IIIIBI 
The Pt distributions for /x"'" and /i^ differ due to the CP- 
breaking valance quark contribution to the proton parton 
distributions. However, given the high CoM energy of 
LHG collisions, this is a small effect. We thus do not show 
the /i+ Pt distributions, being nearly indistinguishable 
from those of /i^. 

In order to get accurate statistics in regions with low 
numbers of events, all of our simulations except those 
of Z + jet have at least four times as many simulated 
events as expected LHG events for the given integrated 
luminosity. For the Z + jet background there are ~ 1.3 
simulated events per LHG event. All event numbers are 
then rescaled for our figures, to the number of events 
expected at the LHG. 

All backgrounds have been significantly reduced by our 
choice of cuts, and the remaining background is dom- 
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inated by Zvv. This can be understood given that it 
passes Z selection criteria, and contains a pair of neu- 
trinos, implying larger amounts of fj^ than other back- 
grounds and a greater resilience to our missing energy 
cut. 

Fig. |9] shows event numbers corresponding to an inte- 
grated luminosity of 5fb~^ and a CoM energy of 7TeV. It 
is clear that our selection cuts are effectively distinguish- 
ing signal from background. Results for niy^ = lOGeV 
and TO^ = 700GeV indicate an excess of a few events af- 
ter integration across all bins, demonstrating the poten- 
tial for constraining the model using current data. For 
= 30GeV the signal strength is significantly weaker. 
In this case, a CoM energy of 8TeV and an integrated 
luminosity of 15fb~^ are required for the signal to be vis- 
ible, as demonstrated in Fig. [TOl With the intention of 
studying heavier dark matter masses, we focus primarily 
on the higher integrated luminosities and CoM energies, 
for which the expected signal is significantly enhanced. 

Looking to the future, we turn our attention to the 
LHC's design CoM energy of 14TeV, and a larger inte- 
grated luminosity of lOOfb^^. Figs. [11] - [12] show the 
expected signal for varying model parameters and 
m^. Signal to background in most bins for the regions 
of parameter space considered exceed 10%, reaching up 
to 70% in some bins for = 30GeV and = ITeV. 
We expect such an excess, if it exists, to be clearly visi- 
ble in future data. Consequently, the non-observation of 
this signal has the potential to place strong constraints 
on this type of model. 

The signal strength decreases as the dark matter mass 
increases, approaching invisibility as approaches 
lOOGeV. Though not implemented here, this decrease in 
signal to background with increasing could poten- 
tially be offset by optimizing cuts on measured events 
for a given dark matter mass. Alternatively, one could 
relax the requirement that x constitute all of the dark 
matter, allowing an increase in which would in turn 
scale the production cross section. An example of such 
non-saturating dark matter can be seen in Fig. 1131 Here, 
]ud is taken to be 3, in contrast to the value of 1.67 re- 
quired to satisfy the relic abundance constraint. With 
jud = 3, the X contributes only ^ 10% to the total dark 
matter abundance. 



ticular, we advance the process — > XX-^; where the Z 
decays muonically. The final state appears as a pair of 
high j)T muons, with an invariant mass in the Z window, 
and large amounts of -^y. Despite having a cross section 
which is low relative to processes with jets in the final 
state, this process has few Standard Model backgrounds, 
and these can be tamed by modest event cuts. By ap- 
plying a specific model for dark matter, with direct dark 
matter coupling to quarks, we have demonstrated that 
this process should be highly visible in future data from 
LHC upgrades. We have found that for certain choices 
of model parameters an excess of events should be visible 
across a broad range of energy bins; some bins contain 
up to a 70% signal to background ratio. As a result, this 
mono-Z search can at the very least provide important 
complementary information to jet and photon based dark 
matter searches. This mono-Z search is relevant whether 
or not its discovery potential is competitive with mono- 
jets and mono-photons, due to the unique kinematical 
aspects of Z decay in distinguishing mono-Z's from other 
bremsstrahlung based search channels. 

ACKNKOWLEDGEMENTS 

We thank Elisabetta Barberio, Bhaskar Dutta, Teruki 
Kamon, Martin White, Antonio Limosani and Nicholas 
Stetzer for helpful discussions. NFB was supported by 
the Australian Research Council, AJG was supported by 
the Commonwealth of Australia, TJW was supported in 
part by U.S. DoE grant DEFG0585ER40226, and LMK, 
JBD, and TDJ are supported in part by the U.S. DoE 
grant DE-SC0008016. LMK acknowledges the hospital- 
ity of Australian National University, and the Research 
School of Astronomy and Astrophysics there, where some 
of this work was carried out. 



VI. CONCLUSIONS 



There are many different search channels for dark mat- 
ter at the LHC, most of those being dependent on the 
specifics of the model under consideration, supersymmet- 
ric or otherwise. A key hallmark of any of these searches 
are some set of high events, against whose momentum 
large amounts of $t can be reconstructed. A complete 
dark matter search must take into account signatures in 
all possible channels. In this work, we have pointed out 
the relevance of mono-Z (electroweak bremsstrahlung) 
processes as a search channel for dark matter. In par- 
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. 9. Events passing selection criteria in Section IIII Bl at — 7TeV and 5/6 ^ of data, for m,, — 700GeV, comparing 
= lOGeV and 30GeV. 
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FIG. 10. As Fig[9l with y/s = 8TeV and 15fb"^ of data, for = 30GeV and = 700GeV. 
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FIG. 11. As Fig.H with ^ = 14TeV and lOOfb"^ of data, for = 30GeV, comparing m„ = 700GeV and lOOOGeV. 
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FIG. 12. As Fig. H with y/s = 14TeV and lOOfb"^ of data, for = 700GeV, comparing =50 GeV, and 100 GeV. 
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FIG. 13. As Fig. 1121 with — lOOGeV and = 700GeV, for fud ~ 3, which corresponds to 10% of present relic abundance. 
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